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Abstract 



We introduce a new approach to computing an approximately maximum s-t flow in a capaci- 
tated, undirected graph. This flow is computed by solving a sequence of electrical flow problems. 
Each electrical flow is given by the solution of a system of linear equations in a Laplacian matrix, 
and thus may be approximately computed in nearly-linear time. 

Using this approach, we develop the fastest known algorithm for computing approximately 
maximum s-t flows. For a graph having n vertices and m edges, our algorithm computes a (1— e)- 
approximately maximum s-t flow in timd^j O (mn 1 /^ -11 / 3 ). A dual version of our approach 

computes a (f + e)-approximately minimum s-t cut in time O (m + n 4 / 3 e~ 8 / 3 ) , which is the 
fastest known algorithm for this problem as well. Previously, the best dependence on m and 
n was achieved by the algorithm of Goldberg and Rao (J. ACM 1998), which can be used to 
compute approximately maximum s-t flows in time O (m^/ne _1 ), and approximately minimum 
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1 We recall that O (/(m)) denotes 0(f(m) log c f(m)) for some constant c. 




1 Introduction 



The maximum s-t flow problem and its dual, the minimum s-t cut problem, are two of the most 
fundamental and extensively studied problems in Operations Research and Optimization |1 6|, [JJ . 
They have many applications (see [2]) and are often used as subroutines in other algorithms (see [HI 
[T7]). Many advances have been made in the development of algorithms for this problem (see 
Goldberg and Rao [TT] for an overview). However, for the basic problem of computing or (1 — e)- 
approximating the maximum flow in undirected, unit-capacity graphs with m = 0(n) edges, the 
asymptotically fastest known algorithm is the one developed in 1975 by Even and Tarjan [3], which 
takes time 0(n 3 / 2 ). Despite 35 years of extensive work on the problem, this bound has not been 
improved. 

In this paper, we introduce a new approach to computing approximately maximum s-t flows 
and minimum s-t cuts in undirected, capacitated graphs. Using it, we present the first algorithms 
that break the 0(n 3//2 ) complexity barrier described above. In addition to being the fastest known 
algorithms for this problem, they are simple to describe and introduce techniques that may be 
applicable to other problems. In them, we reduce the problem of computing maximum flows 
subject to capacity constraints to the problem of computing electrical flows in resistor networks. 
An approximate solution to each electrical flow problem can be found in time O (m) using recently 
developed algorithms for solving systems of linear equations in Laplacian matrices |13| ITS] . 

There is a simple physical intuition that underlies our approach, which we describe here in the 
case of a graph with unit edge capacities. We begin by thinking of each edge of the input graph as 
a resistor with resistance one, and we compute the electrical flow that results when we send current 
from the source s to the sink t. These currents obey the flow conservation constraints, but they 
may not respect the capacities of the edges. To remedy this, we increase the resistance of each edge 
in proportion to the amount of current flowing through it — thereby penalizing edges that violate 
their capacities — and compute the electrical flow with these new resistances. 

After repeating this operation O (m 1 / 3 • poly(l/e)) times, we will be able to obtain a (1 — e)- 
approximately maximum s-t flow by taking a certain average of the electrical flows that we have 
computed, and we will be able to extract a (1 + e)-approximately minimum s-t cut from the vertex 
potentials 2 ! This will give us algorithms for both problems that run in time O (m 4 / 3 • poly(l/e)). 
By combining this with the graph smoothing and sampling techniques of Karger [p2] , we can get 
a (1 — e)-approximately maximum s-t flow in time O (mn'/V 11 ' 3 ). Furthermore, by applying 
the cut algorithm to a sparsifier |4. of the input graph, we can compute a (1 + e)-approximately 
minimum s-t cut in time O (m + n 4 / 3 e~ 8 / 3 ). 

We remark that the results in this paper immediately improve the running time of algorithms 
that use the computation of an approximately maximum s-t flow on an undirected, capacitated 
graph as a subroutine. For example, combining our work with that of Sherman [T7] allows us to 
achieve the best currently known approximation ratio of 0{yJ\ogn) for the sparsest cut problem in 
time O (m + n 4 / 3 ) . 

We are hopeful that our approach can be extended to directed graphs and can also eventually 
lead to an algorithm that approximately solves the maximum flow problem in nearly-linear time. 



For clarity, we will analyze these two cases separately, and they will use slightly different rules for updating the 
resistances. 
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1.1 Previous Work on Maximum Flows and Minimum Cuts 

The best previously known algorithms for the problems studied here are due to Goldberg and Rao. 
In a breakthrough paper, Goldberg and Rao [TT] developed an algorithm for computing exact maxi- 
mum s-t flows in directed or undirected capacitated graphs in time 0(m min(n 2//3 , m 1 / 2 ) log(n 2 /m) log U), 
assuming that the edge capacities are integers between 1 and U. When we are interested in find- 
ing (1 — e)-approximately maximum s-t flow, the dependence on log U can be removed and, by 
employing the smoothing technique of Karger [12] , one can obtain a running time of 

O (rn^Jne~ l ) . 

By applying their algorithm to a sparsifier, as constructed by Benczur and Karger [3], Goldberg 
and Rao show how to compute a (1 + e)-approximately minimum s-t cut in an undirected graph in 
time 

O (m + n 3/ V 3 ) . 

Their work was the culmination of a long line of papers on the problem; we refer the reader to 
their paper for an extensive survey of earlier developments in algorithms for computing maximum 
s-t flows. In more recent work, Daitch and Spielman [7] showed that fast solvers for Laplacian 
linear systems |1 8|, [T3] could be used to make interior-point algorithms for the maximum flow 
and minimum-cost flow problems run in time O (m 3//2 log J7) , and Madry [T3j showed that one 
can approximate a wide range of cut problems, including the minimum s-t cut problem, within a 
polylogarithmic factor in almost linear time. 



1.2 Outline 

We begin the technical part of this paper in Section [2] with a review of maximum flows and electrical 
flows, along with several theorems about them that we will need in the sequel. In Section [3] 
we give a simplified version of our approximate maximum-flow algorithm that has running time 
O (m 3//2 e~ 5 / 2 ). In Section |4j we will show how to improve the running time of our algorithm 
to O (m 4 / 3 e~ 3 ); we will then describe how to combine this with existing graph smoothing and 

) 



1/3,-11/3 



ran ' e 



sparsification techniques to compute approximately maximum s-t flows in time O 
and to approximate the value of such flows in time O (ra + n 4 / 3 e~ 8 / 3 ). In Section [5 we present a 
variant of our algorithm that computes approximately minimum s-t cuts in time O (m + n 4//3 e -8 / 3 ) 



2 Maximum Flows, Electrical Flows, and Laplacian Systems 
2.1 Graph Theory Definitions 

Throughout the rest of the paper, let G = (V, E) be an undirected graph with n vertices and m 
edges. We distinguish two vertices, a source vertex s and a sink vertex t. We assign each edge e a 
nonzero integral capacity u e G Z + , and we let U := max e u e j min e u e be the ratio of the largest to 
the smallest capacities. 

We arbitrarily orient each edge in E; this divides the edges incident to a vertex v GV into the 
set E + (v) of edges oriented towards v and the set E~(v) of edges oriented away from v. These 
orientations are merely for notational convenience. We use them to interpret the meaning of a 
positive flow on an edge. If an edge has positive flow and is in E + (v), then the flow is towards v. 
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Conversely, if it has negative flow then the flow is away from v. One should keep in mind that our 
graphs are undirected and that the flow on an edge can go in either direction, regardless of this 
edge's orientation. 

We now define our primary objects of study, s-t cuts and s-t flows. 

Definition 2.1 (Cuts). An s-t cut is a partition (S, V\S) of the vertices into two disjoint sets such 
that s 6 S and t E V\S. The capacity u(S) of the cut is defined to be the sum u(S) := YleeE(S) 
where E(S) C E is the set of edges with one endpoint in S and one endpoint in V \ S. 

Definition 2.2 (Flows). An s-t flow is a function f : E — )• 1R that obeys the flow-conservation 
constraints 

E E /( e )=° for alive V\{s,t}. 

e€E~(v) e£E+(v) 

The value |/| of the flow is defined to be the net flow out of the source vertex, \f \ := YleeE-(s) f( e )~ 

It follows easily from the flow conservation constraints that the net flow out of s is equal to the 
net flow into t, so |/| may be interpreted as the amount of flow that is sent from s to t. 

2.2 Maximum Flows and Minimum Cuts 

We say that an s-t flow / is feasible if |/(e)| < u e for each edge e, i.e., if the amount of flow 
routed through any edge does not exceed its capacity. The maximum s-t flow problem is that of 
finding a feasible s-t flow in G of maximum value. We denote a maximum flow in G (with the given 
capacities) by /*, and we denote its value by F* := |/*|. We say that / is a (1 — e)-approximately 
maximum flow if it is a feasible s-t flow of value at least (1 — e)F*. 

To simplify the exposition, we will take e to be a constant independent of m throughout the 
paper, and m will be assumed to be larger than some fixed constant. However, our analysis will 
go through unchanged as long as e > f2(m _1//3 ). In particular, our analysis will apply to all e for 
which our given bounds are faster than the 0(m 3 / 2 ) time required by existing exact algorithms. 

One can easily reduce the problem of finding a (1 — e)-approximation to the maximum flow in 
an arbitrary undirected graph to that of finding a (1 — e/2)-approximation in a graph in which 
the ratio of the largest to smallest capacities is polynomially bounded. To do this, one should 
first compute a crude approximation of the maximum flow in the original graph. For example, one 
can compute the s-t path of maximum bottleneck in time 0(m + ralogn) |16| Section 8.6e], where 
we recall that the bottleneck of a path is the minimum capacity of an edge on that path. If this 
maximum bottleneck of an s-t path is B, then the maximum flow lies between B and mB. This 
means that there is a maximum flow in which each edge flows at most mB, so all capacities can be 
decreased to be at most mB. On the other hand, if one removes all the edges with capacities less 
eB /2m, the maximum flow can decrease by at most eB/2. So, we can assume that the minimum 
capacity is at least eB/2m and the maximum is at most Bm, for a ratio of 2m 2 /e. Thus, by a 
simple scaling, we can assume that all edge capacities are integers between 1 and 2m 2 /e. 

The minimum s-t cut problem is that of finding the s-t cut of minimum capacity. The Max 
Flow-Min Cut Theorem ([10, 8j states that the capacity of the minimum s-t cut is equal to F* , 
the value of the maximum s-t flow. 
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In particular, the Max Flow-Min Cut Theorem implies that one can use the capacity of any s-t 
cut as an upper bound on the value of any feasible s-t flow, and that the task of finding the value 
of the maximum flow is equivalent to the task of finding the capacity of a minimal s-t cut. 

One should note, however, that the above equivalence applies only to the values of the flow 
and the capacity and that although one can easily obtain a minimum s-t cut of a graph given its 
maximum flow, there is no known procedure that obtains a maximum flow from minimum s-t cut 
more efficiently than by just computing the maximum flow from a scratch. 



2.3 Electrical Flows and the Nearly Linear Time Laplacian Solver 

In this section, we review some basic facts about electrical flows in networks of resistors and 
present a theorem that allows us to quickly approximate these flows. For an in-depth treatment of 
the background material, we refer the reader to jB]. 

We begin by assigning a resistance r e > to each edge e G E, and we collect these resistances 
into a vector r G IR m . For a given s-t flow /, we define its energy (with respect to r) as 

SrU) :=5> e / 2 (e). 

e 

The electrical flow of value F (with respect to r, from s to t) is the flow that minimizes £ r {f) 
among all s-t flows / of value F. This flow is easily shown to be unique, and we note that it need 
not respect the capacity constraints. 

From a physical point of view, the electrical flow of value one corresponds to the current that 
is induced in G if we view it as an electrical circuit in which each edge e has resistance of r e , and 
we send one unit of current from s to t, say by attaching s to a current source and t to ground. 



2.3.1 Electrical Flows and Linear Systems 

While finding the maximum s-t flow corresponds to solving a linear program, we can compute 
the electrical flow by solving a system of linear equations. To do so, we introduce the edge-vertex 
incidence matrix B, which is an n x m matrix with rows indexed by vertices and columns indexed 
by edges, such that 

r l if eG £?-(«), 
B r , ~- < -1 if e G E+(v), 
otherwise. 

If we treat our flow / as a vector / G IR m , where we use the orientations of the edges to 
determine the signs of the coordinates, the v th entry of the vector B T f will be the difference 
between the flow out of and the flow into vertex v. As such, the constraints that one unit of flow 
is sent from s to t and that flow is conserved at all other vertices can be written as 

B T f = Xs,t, 

where Xs,t is the vector with a 1 in the coordinate corresponding to s, a —1 in the coordinate 
corresponding to t, and all other coordinates equal to 0. 

We define the (weighted) Laplacian L of G (with respect to the resistances r) to be the n x n 
matrix 

L := B CB T , 
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where C is the mxm diagonal matrix with C eje = c e = l/r e . One can easily check that its entries 
are given by 

YueeE+(u)YjE~{u) °e '^U = V, 

-c e if e = (u, v) is an edge of G, and 

1 otherwise. 

Let R = C _1 be the diagonal matrix with R ee = r e . The energy of a flow / is given by 

S r (f):=Y j ref(ef=fRf= R^f ' 



The electrical flow of value 1 thus corresponds to the vector / that minimizes Lr 1/2 / subject to 

Bf = Xs,t- If / is an electrical flow, it is well known that it is a potential flow, which means that 
there is a vector (j> £ TR V such that 

f(u,v) = . 

fu,v 

That is, 

/ = CB T cj> = R l B T cj). 

Applying Bf = Xs,t, we have Bf = BCB T <p = Xs,t, and hence the vertex potentials are given 
by 

<f> = L ] Xs,t, 

where L* denotes the Moore- Penrose pseudo-inverse of L. Thus, the electrical flow / is given by 
the expression 

/ = CB T tf X s,t- 
This lets us rewrite the energy of the electrical flow of value 1 as 

£ r (f)=f T Rf = (xs/ti T BC T ) R (cB T P X s,t) = Xs,t^LL^ X s,t = Xs/^ X s,t = <t> T Lct>. 

(1) 



2.3.2 Effective s-t Resistance and Effective s-t Conductance 

Our analysis will make repeated use of two basic quantities from the theory of electrical networks, 
the effective s-t resistance and effective s-t conductance. 

Let / be the electrical s-t flow of value 1, and let (f> be its vector of vertex potentials. The 
effective s-t resistance of G with respect to the resistances r is given by 

R cS (r) = <f>(s) - 4>(t). 

Throughout this paper, we will only look at the effective resistance between the vertices s and t in 
the graph G, so we suppress these letters in our notation and simply write R e f[(r). 

Using our linear algebraic description of the electrical flow and Equation ([!]) , we have 

Res(r) = 4>(s) - Ht) = Xs,t T <f> = Xs,t T tf X s,t = £r(f ). 

This gives us an alternative description of the effective s-t resistance as the energy of the electrical 
flow of value 1 . 
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It will sometimes be convenient to use the related notion of the effective s-t conductance of G 
with respect to the resistances r, which we define by 

C cS (r) = l/R eS (r). 

We note that this is equal to the value of the electrical flow in which (j)(s) — <j>{t) = 1. 



2.3.3 Approximately Computing Electrical Flows 

From the algorithmic point of view, the crucial property of the Laplacian L is that it is symmetric 
and diagonally dominant, i.e., for any u, Y2v'^u — ^u,v This allows us to use the result of 

Koutis, Miller, and Peng [13], which builds on the work of Spielman and Teng |18] . to approximately 
solve our linear system in nearly-linear time. By rounding the approximate solution to a flow, we 
can prove the following theorem (see Appendix |A| for a proof). 

Theorem 2.3 (Fast Approximation of Electrical Flows). For any 5 > 0, any F > 0, and any 

vector r of resistances in which the ratio of the largest to the smallest resistance is at most R, we 
can compute, in time O (mlogR/S), a vector of vertex potentials cj) and an s-t flow f of value F 
such that 

a. £ r {f) < (1 + S)£ r (f), where f is the electrical s-t flow of value F, and 

b. for every edge e, 

where f is the true electrical flow. 

c. 

We will refer to a flow meeting the above conditions as a 5-approximate electrical flow. 



2.4 How the Resistance of an Edge Influences the Effective Resistance 

In this section, we will study how changing the resistance of an edge affects the effective resistance 
of the graph. This will be a key component of the analysis of our O (m 4 / 3 • poly(l/e)) algorithms. 
We will make use of the following standard fact about effective conductance; we refer the reader to 
j6l Chapter IX. 2, Corollary 5] for a proof. 

Fact 2.4. For any G = (V, E) and any vector of resistances r, 

C c{i (r) = mm > . 

0t=o (u,v)eE (u ' v > 

Furthermore, the equality is attained for <f> being vector of vertex potentials corresponding to the 
electrical s-t flow of G (with respect to r) of value 1/R e s(r). 

Corollary 2.5 (Rayleigh Monotonicity) . If r' e > r e for all e£fi, then R c g(r r ) > R e s(r). 
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Proof. For any 0, 

0» ~ 0^) 2 < ^ (<Pu ~ 4>v) 2 

(u,v)eE r {u,v) ( u ,v)eE r ( u > v ) 

so the minima of these expressions over possible values of eft obey the same relation, and thus 
Ces(f') < C e fj(r). Inverting both sides of this inequality yields the desired result. □ 

Our analysis of the O (m 4 / 3 ) algorithm will require the following lemma, which gives a lower 
bound on the effect that increasing the resistance of an edge can have on the effective resistance. 

Lemma 2.6. Let f be an electrical s-t flow on a graph G with resistances r. Suppose that some 
edge h = (i,j) accounts for a (3 fraction of the total energy of f , i.e., 

f(h) 2 r h = (3£ r (f). 

For some 7 > 0, define new resistances r' such that r' h = 777^ and r' e = r e for all e//i. Then 

2 ^.fl *«<->■ 

In particular: 

• If we "cut" the edge h by setting 7 = 00, then 

R c s(r) 



Res(r') > 



1-/3 



• // we slightly increase the effective resistance of the edge h by setting 7 = (1 + e) with e < 1, 
then 

Res(r') > - n — , 1+ ^ 7wR c s(r) > ( 1 + — ^ R cS (r). 

Proof. The assumptions of the theorem are unchanged if we multiply / by a constant, so we may 
assume without loss of generality that / is the electrical s-t flow of value 1/R c s(r). If <f>' is the 
vector of vertex potentials corresponding to /, this gives <f>{ — <j){ = 1. Since adding a constant to 



the potentials doesn't change the flow, we may assume that (f> s = 1 and 4>t = 0. By Fact 2.4 



n , , sr (4-d) 2 , x - {4-d) 2 

Ceff(r) = 2^ = + ■ 

(u,v)EE r ^ Th (u,v)GE\{h} r ^ 

The assumption that h contributes a /3 fraction of the total energy implies that, in the above 
expression, 



if _ 2 



and thus 



/3C eff (r), 



E — r J - = (l-P)C eS (r). 

{u,v)eE\{h} ( u>v ) 
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We will obtain our bound on C e g(r') by plugging the original vector of potentials $ into the 



expression in Fact 2.4 



C eS (r = mm } '— < } '— 

4>t=o (u,v)eE («> u ) (u,v)eE 

_ {4 - rt? | ^ (d-ti) 2 = , y (^-^) 2 

h {u,v)dE\{h} ( u > v ) ' n {u,v)£E\{h} (u > v) 

- C eS (r) + (l-(3)C cS (r) = C c(X (r' 1 ' ° 
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Since R e ff(r) = 1/C c g(r) and R e g(r') = l/C e g(r'), the desired result follows. □ 

3 A Simple O (m 3/ \ 5//2 )-Time Flow Algorithm 

Before describing our O (m 4 / 3 e -3 ) algorithm, we will describe a simpler algorithm that finds a 
(1 — e)-approximately maximum flow in time O (m 3 / 2 e~ 5 / 2 ). Our final algorithm will be obtained 
by carefully modifying the one described here. 

The algorithm will employ the multiplicative weights update method [3j 115] . In our setting, 
one can understand the multiplicative weights method as a way of taking an algorithm that solves 
a flow problem very crudely and, by calling it repeatedly, converts it into an algorithm that gives 
a good approximation for the maximum flow in G. The crude algorithm is called as a black-box, 
so it can be thought of as an oracle that answers a certain type of query. 

In this section, we provide a self-contained description of the multiplicative weights method 
when it is specialized to our setting. In Section [3. 1[ we will describe the requirements on the oracle, 
give an algorithm that iteratively uses it to obtain a (1 — e)-approximately maximum flow, and 
state how the number of iterations required by the algorithm depends on the properties of the 



oracle. In Section 3.2 we will describe how to implement the oracle using electrical flows. Finally, 



in Section 3.3 we will provide a simple proof of the convergence bound set forth in Section [3TT 



3.1 Multiplicative Weights Method: From Electrical Flows to Maximum Flows 

For an s-t flow /, we define the congestion of an edge e to be the ratio 

, v |/(e)| 
cong / (e) := 

u e 

between the flow on an edge and its capacity. In particular, an s-t flow is feasible if and only if 
congj(e) < 1 for all e £ E. 

The multiplicative weights method will use a subroutine that we will refer to as an (e, p)-oracle. 
This oracle will take as input a number F and a vector w of edge weights. For any F < F* , 
we know that there exists a way to route F units of flow in G so that all of the edge capacities 
are respected. Our oracle will provide a weaker guarantee: When F < F*, it will satisfy all of 
the capacity constraints up to a multiplicative factor of and it will satisfy the average of these 



J Up to polynomial factors in 1/e, the value of p will be Q(y/m) in this section, and 0(m 1,/3 ) later in the paper 
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constraints, weighted by the Wi, up to a (much better) multiplicative factor of (1 + e). When 
F > F* , the oracle will either output an s-t flow satisfing the conditions above, or it will return 
"fail" . 

Formally, we will use the following definition: 

Definition 3.1 ((e, p) oracle). For e > and p > 0, an (e,p) oracle is an algorithm that, given a 
real number F > and a vector w of edge weights with w e > 1 for all e, returns an s-t flow f such 
that: 

1- If F < F* , then it outputs an s-t flow f satisfying: 

(i) l/l = F; 

(ii) ^e^econgy^e) < (1 + e)|u;|i, where \ w\i := ]P g w e ; 
(Hi) max e congj(e) < p. 

2. If F > F* , then it either outputs a flow f satisfying conditions (i), (ii), (Hi) or outputs 
"fail". 

Our algorithm will be given a flow value F as an input. If F < F* , it will return a flow of value 
at least (1 - 0(e))F. If F > F*, it will either return a flow of value at least (1 - 0(e)) F (which 
may occur if F is only slightly greater than F*) or it will return "fail". This allows us to find a 



(1 — 0(e))-approximation of F* using binary search. As outlined in Section 2.2 we can obtain a 
crude bound B in time 0(m + nlogn) such that B < F* < mB, so the binary search will only call 
our algorithm a logarithmic number of times. 

In Figure [T] we present our simple algorithm, which applies the multiplicative weights update 
routine to approximate the maximum flow by calling an (e, p)-flow oracle. The algorithm initializes 
all of the weights to 1 and calls the oracle with these weights. The call returns a flow that satisfies 
conditions (i), (ii), and (iii) defined above. It then multiplies the weight of each edge e by (1 + 
^congji(e)). Note that if the congestion of an edge is poor, say close to p, then its weight will 
increase by a factor of (1 + e). On the other hand, if the flow on an edge is no more than its 
capacity, then the new weight of the edge is essentially unchanged. This will put a larger fraction 
of the weight on the violated constraints, so the oracle will be forced to return a solution that 
comes closer to satisfying them (possibly at the expense of other edges). In the end, we return the 
average of all of the flows as our answer. 

The key point in analyzing this algorithm is that the total weight on G does not grow too 
quickly, due to the average congestion constraint on the flows returned by the oracle O. However, 
if an edge e consistently suffers large congestion in a sequence of flows returned by O, then its 
weight increases rapidly relative to the total weight, which will significantly penalize any further 
congestion the edge in the subsequent flows. If this were to occur too many times, its weight would 



exceed the total weight, which obviously cannot occur. In Section 3.3, we will prove the following 
theorem by showing that our algorithm converges in 2/olnm/e 2 iterations. 

Theorem 3.2 (Approximating Maximum Flows by Multiplicative Weights). For any < e < 1/2 

and p > 0, given an (e,p)-flow oracle with running time T(m, 1/e, U), one can obtain an algorithm 
that computes a (1 — 0(e))- approximate maximum flow in a capacitated, undirected graph in time 
d(pe- 2 -T(m,l/e,U)). 
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Input : A graph G = (V,E) with capacities {u e } e , a target flow value F, and an (e, p)-oracle O 
Output: Either a flow J, or "fail" indicating that F > F*; 

Initialize w° e <- 1 for all edges e, and N 4- 
for i := 1, . . . , N do 

Query with edge weights given by w % ~ 1 and target flow value F 

if returns "fail" then return "fail" 

else 

Let f l be the returned flow 

w\ <— k4 -1 (1 + ^cong^i(e)) for each e <E E 
end 
end 

return / <-^(Ei/') 

Figure 1: Multiplicative-weights-update routine 



Note that the number of iterations of the algorithm above grows linearly with the value of p, 
which we call the width of the oracle. Intuitively, this should be necessary because the final flow 
across an edge e is equal to the average of the flows sent over it by all /* . If we send p ■ u e units of 
flow across the edge e in some step, then we will need at least iterations to drop the average 
to 10 

3.2 Constructing an Oracle of Width 3-^/m/e Using Electrical Flows 



Given Theorem 3.2 our problem is thus reduced to designing an efficient oracle that has a small 

m log e^ 1 ) time implementation of an (e, 3y/m/e)- 



3.2 



width. In this subsection, we will give simple O 
flow oracle for any < e < 1/2. By Theorem 
time algorithm for finding an approximately maximum flow. 
To build such an oracle, we set 



this will immediately yield an O {rrftl^e 5//2 ) 



1 f : e\w\i 
3m 



2 i ™e + -br ) ( 2 ) 



for each edge e, and we use the procedure from Theorem 2.3 to approximate the electrical flow 
that sends F units of flow from s to t in a network whose resistances are given by the r e . The 
pseudocode for this oracle is shown in Figure [2} 



We now show that the resulting flow / has the properties required by Definition 3.1 Since 
| /| = F by construction, we only need to demonstrate the bounds on the average congestion 
(weighted by the w e ) and the maximum congestion. We will use the basic fact that electrical flows 
minimize the energy of the flow. Our analysis will then compare the energy of / with that of an 
optimal max flow. Intuitively, the w e term in Equation Q guarantees the bound on the average 
congestion, while the e|u>|i/(3m) term guarantees the bound on the maximum congestion. 

4 Strictly speaking, it is possible that we could do better than this by sending a large amount of flow across the 
edge in the opposite direction. However, nothing in our algorithm aims to obtain this kind of cancellation, so we 
shouldn't expect to be able to systematically exploit it. 
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Input : A graph G = (V,J5) with capacities {u e } et a target flow value F, and edge weights {w e } e 


niifnnt- Fi-t-hoK- :a flr»\n/ f r»r "fail" i nrJ i<~^-|-i n rr F Z7* 

i^/uxjjux . citner a now j , or Tail inaicanng inat r s> r 






re^+z (w e + ^) for each e e £ 






Find an (e/3)-approximate electrical flow / using Theorem 


2.3 


on G with resistances r and target 


flow value F 






if £ r (f) > (1 + e)|«?|i then return "/aii" 






else return / 






Figure 2: A simple implementation of an fe, 3\/m/e\ oracle 



Suppose /* is a maximum flow. By its feasibility, congy*(e) < 1 for all e, so 



< 



1 + Hi. 



Since the electrical flow minimizes the energy, £ r (f*) is an upper bound on the energy of the 



electrical flow of value F whenever F < F*. In this case, Theorem 2.3 implies that our (e/3)- 
approximate electrical flow satisfies 

£r(f) < (i + 1) £r(f*) < (i + 1) 2 Hi < (1 + e ) Hi- (3) 

This shows that our oracle will never output "fail" when F < F* . It thus suffices to show that 
the energy bound £ r (f) > (l + e )Hii which holds whenever the algorithm does not return "fail", 
implies the required bounds on the average and worst-case congestion. To see this, we note that 
the energy bound implies that 

Y,w e (congj(e)) 2 < (l + e)Hi, (4) 

e 

and, for all e € E, 



e\w\i 
3m 

By the Cauchy-Schwarz inequality, 



congj(e)) <(l + e)Hi. (5) 



^2 w e congj(e)^ < |tw|i w e (congj(e) 



(6) 
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so Equation Q gives us that 

^ w e congjr(e) < Vl + e|w|i < (1 + e) | w|i, (7) 

e 

which is the required bound on the average congestion. Furthermore, Equation ^ and the fact 
that e < 1/2 implies that 



3m(l + e) 



congj(e) < y ; < 3yW 



for all e, which establishes the required bound on the maximum congestion. So our algorithm 
implements an (e, 3-y/m7e)-oracle, as desired. 

To bound the running time of this oracle, recall that we can assume all edge capacities lie 
between 1 and U = m 2 /e and compute 

„ r e T r,max e r e + e \w\i , T9 m + e „ /. , , n ni\ , x 

J R = max— <C/ 2 f , 1 11 < U 2 = O I m/e . (8 

e,e' r e / e e \ J 

This establishes an upper bound on the ratio of the largest resistance to the smallest resistance. 
Thus, by Theorem 2.3 the running time of this implementation is O (mlogR/e) = O (mlogl/e). 
Combining this with Theorem 3.2 we have shown 

Theorem 3.3. For any < e < 1/2, the maximum flow problem can be (1 — e)- approximated in 
O (m 3 / 2 e~ 5 / 2 ) time. 

3.3 The Convergence of Multiplicative Weights 

In this section, we prove Theorem |3.2| by analyzing the multiplicative weights update algorithm 
shown in Figure [TJ 

Our analysis will be based on the potential function /ij := |tu*|i. Clearly, [xq = m and this 
potential only increases during the course of the algorithm. It follows from condition (l.ii) of 
Definition 3.1 that if we run the (e, p)-oracle O with F = F*, then 

5^u4cong /i+ i(e) < (1 +e)|to i |i, for all i > 1. (9) 



By condition (l.iii) of Definition 3.1 



congji(e) < p, for alH > 1 and any edge e. (10) 
We start by upper bounding the total growth of fa thoughout the algorithm. 
Lemma 3.4. For any i > 0, 

< Mi ex P 



p 

, N \U = „.,r < mP^n ( ^^ n) =n°( 1 / £ ). 



In particular, \\w ||i = < mexp ( - 
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Proof. For any i > 0, we have 



Pi 

e e x ' ' e ' & 

where the last inequality follows from Q. Thus, we can conclude that 

(l + e)e, ,. , (l + e)e, fO- + e)e 

IH+l < Pi H to i = Pi(l H ) < Piexp 

P P VP 

as desired. The lemma follows. □ 

One of the consequences of the above lemma is that whenever we make a call to the or- 
acle, the total weight is at most n 0<yl ^\ Thus, the running time of our algorithm is 
O (pe~ 2 • T(m,l/e, U,m 0{ - l/e) )) as claimed. 

Next, we bound the final weight of a particular edge e with the congestion congj(e) that 

this edge suffers in our final flow /. 
Lemma 3.5. For any edge e and i > 0, 

>e > exp ( — — — ^congp (e) ] . 



IV, 



p m 



In particular, > exp ^ ^\i- t )p con Ej( e ) 
Proof. For any i > 0, we have 



f[(l+ ^cong /3 (e)) > [J exp (ii-^ C ong p (e)) , 



j>1 v ' " j>i 



where we used (10) and that for any 1/2 > e > and x G [0, 1]: 

exp((l — e)ex) < (1 — ex). 

Now, the lemma follows since 

>jj /(l-e)e \ /(i_ e ) e * \ 
«e ^ 1 1 exp I congp (e) = exp 2^ con gp( e ) > 

and for i = N 

( Mg y , A /(l + e)eiV A 

w e > exp I — - — 2^ cor| g/n e ) I = ex P ^ (i _ e ) p con g7v e ) J ■ 



□ 
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Finally, by Lemmas |3.4| and 3.5, we conclude that for any edge e 
(l + e)e7V 



m exp 
This implies that 



> 



PN 



\w \i > w e > exp 



1 + e)eN 



cong 7 (e) 



(1 — e)plnm 



e(l-e) 



for every edge e. Thus, we see that / is a feasible s-t flow and, since each p has throughput i 7 * 
the throughput |/| of f is j^-F* > (1 - 0(e))F* for 1/2 > e > 0, as desired. 



4 An O (mn 1//3 e 1X / 3 ) Algorithm for Approximate Maximum Flow 

In this section, we modify our algorithm to run in time O (m 4//3 e~ 3 ) . We then combine this with the 
smoothing and sampling techniques of Karger [12] to obtain an O (mn 1 / 3 e _11 / 3 )-time algorithm. 

For fixed e, the algorithm in the previous section required us to compute O (m 1 / 2 ) electrical 
flows, each of which took time O (m), which led to a running time of O (m 3 / 2 ). To reduce this to 
we'll show how to find an approximate flow while computing only O (m 1 / 3 ) electrical 

flows. 



Our analysis of the oracle from Section 3.2 was fairly simplistic, and one might hope to improve 
the running time of the algorithm by proving a tighter bound on the width. Unfortunately, the 
graph in Figure [3] shows that our analysis was essentially tight. The graph consists of k parallel 
paths of length k connecting s to t, along with a single edge e that directly connects s to t. The 
max flow in this graph is k + 1. In the first call made to the oracle by the multiplicative weights 
routine, all of the edges will have the same resistance. In this case, the electrical flow of value k + 1 
will send (k + l)/2k units of flow along each of the k paths and (k + l)/2 units of flow across e. 
Since the graph has m = Q(k 2 ), the width of the oracle in this case is 0(m 1//2 ). 



4.1 The Improved Algorithm 

The above example shows that it is possible for the electrical flow returned by the oracle to exceed 
the edge capacities by 0(m 1//2 ). However, we note that if one removes the edge e from the graph 
in Figure [3] the electrical flow on the resulting graph is much better behaved, but the value of the 
maximum flow is only very slightly reduced. This demonstrates a phenomenon that will be central 
to our improved algorithm: while instances in which the electrical flow sends a huge amount of 
flow over some edges exist, they are somewhat fragile, and they are often drastically improved by 
removing the bad edges. 

This motivates us to modify our algorithm as follows. We'll set p to be some value smaller than 
the actual worst-case bound of O (m 1//2 ). (It will end up being O (m 1//3 ).) The oracle will begin 
by computing an electrical flow as before. However, when this electrical flow exceeds the capacity 
of some edge e by a factor greater than p, we'll remove e from the graph and try again, keeping all 
of the other weights the same. We'll repeat this process until we obtain a flow in which all edges 
flow at most a factor of p times their capacity (or some failure condition is reached), and we'll use 
this flow in our multiplicative weights routine. When the oracle removes an edge, it is added to a 
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k paths 
of length k 



1 edge 
e 




t 



Figure 3: A graph on which the electrical flow sends approximately y/rn units of flow across an 
edge when sending the maximum flow F* from s to t. 



set H of forbidden edges. These edges will be permanently removed from the graph, i.e., they will 
not be included in the graphs supplied to future invocations of the oracle. 

In Figures [4] and [5] we present the modified versions of the oracle and overall algorithm, where 
we have highlighted the parts that have changed from the simpler version shown in Figures [T] and [2] 

4.2 Analysis of the New Algorithm 

Before proceeding to a formal analysis of the new algorithm, it will be helpful to examine what is 
already guaranteed by the analysis from Section [3j and what we'll need to show to demonstrate 
the algorithm's correctness and bound its running time. 

We first note that, by construction, the congestion of any edge in the flow / returned by 
the modified oracle from Figure [4] will be bounded by p. Furt hermore, it enforces the bound 



£r{f) < (1 + e)|to|i; by Equations Q, (JgJ) , and ([7]) in Section 3.2, this guarantees that / will meet 
the weighted average congestion bound required for a (e, p)-oracle. So, as long as the modified 
oracle always successfully returns a flow, it will function as an (e, /o)-oracle, and our analysis from 
Section [3] will show that the multiplicative update scheme employed by our algorithm will yield an 
approximate maximum flow after O (p) iterations. 

Our problem is thus reduced to understanding the behavior of the modified oracle. To prove 
correctness, we will need to show that whenever the modified oracle is called with F < F*, it will 
return some flow / (as opposed to returning "fail"). To bound the running time, we will need to 
provide an upper bound on the total number of electrical flows computed by the modified oracle 
throughout the execution of the algorithm. 

To this end, we will show the following upper bound on the cardinality \H\ and the capacity 
u(H) of the set of forbidden edges, whose proof we postpone until the next section: 
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Input : A graph G = {V,E) with capacities {u e } e , a target flow value F, edge weights {u> e } £ 

and a set H of forbidden edges 
Output: Either a flow / and a set H of forbidden edges, or "fail" indicating that F > F* 

m 

re^+i (w e + 4^) for each e e E \H 
Find an approximate electrical flow / using Theorem 

target flow value F, and parameter S = e/3. 
if £ r (f) > (1 + e)|f|i or s and t are disconnected in Gh then return "fail" 
if there exists e with congj(e) > p then add e to H and start over 

return / 

Figure 4: The modified oracle O' used by our improved algorithm 



2.3 



on := (V,E\H) with resistances r, 



Input : A graph G = (V,E) with capacities {u e } e , and a target flow value F; 
Output: Either a flow J, or "fail" indicating that F > F*\ 

Initialize io° <- 1 for all edges e, H % , 8" 1/3 in 1/3 "» _ an d TV ^ 

for i := 1, . . . , N do 

Query 0' with edge weights given by w % ~ x , target flow value F, and forbidden edge set H 

if returns "fail" then return "fail" 
else 

Let f l be the returned answer 

Replace H with the returned (augmented) set of forbidden edges 

w l e uig -1 (l + 1 cong fi (e)) for each e <E E 
end 
end 

return J 

Figure 5: An improved (1 — 0(e))-approximation algorithm for the maximum flow problem 
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Lemma 4.1. Throughout the execution of the algorithm, 



30m In m 

H < 2 2 

e A p A 

and 

, , 30mFlnm 
< H ) < e 2 p 3 ■ 



If we plug in the value p = (8m 1//3 In 1 / 3 m)/e used by the algorithm, Lemma 4.1 gives the 
bounds \H\ < g(mliim) 1 / 3 and u(iT) < ^eF < eF/12. 

Given the above lemma, it is now straightforward to show the following theorem, which estab- 
lishes the correctness and bounds the running time of our algorithm. 

Theorem 4.2. For any < e < 1/2, if F < F* the algorithm in Figure^will return a feasible s-t 
flow f of value \ f\ = (1 — 0(e)) F in time O (m 4 / 3 e~ 3 ). 

Proof. To bound the running time, we note that, whenever we invoke the algorithm from Theorem 



2.3 , we either advance the number of iterations or we increase the cardinality of H, so the number 
of linear systems we solve is at most N + \H\ < N + ^|(mlnm) 1//3 . 



Equation ([8| implies that the value of R from Theorem 2.3 is 0((m/e)°^), so solving each 



linear system takes time at most O (mlog 1/e). This gives an overall running time of 

0((N+ if(mlnm) 1 / 3 ) m) = O (m^V 3 ' 



as claimed. 

It thus remains to prove correctness. For this, we need to show that if F < F*, then the oracle 
does not return "fail", which would occur if we disconnect s from t or if S r (f) > (1 + e)|u>|i. By 



Lemma 4.1 and the comment following it, we know that throughout the whole algorithm Gh has 
maximum flow value of at least F* — eF/12 > (1 — e/12) F and thus, in particular, we will never 
disconnect s from t. 

Furthermore, this implies that there exists a feasible flow in our graph of value (1 — e/12)F, 
even after we have removed the edges in H. There is thus a flow of value F in which every edge has 
congestion at most 1/ (1 — e/12). We can therefore use the argument from Section [3.2| (Equation ^ 
and the lines directly preceding it) to show that we always have 

£ r (J) < (1 + e/12) 2 (l + e/3) 2 |w|i < (1 + 
as required. □ 



The above theorem allows us to apply the binary search strategy that we used in Section 3.1 
This yields our main theorem: 

Theorem 4.3. For any < e < 1/2, the maximum flow problem can be (1 — e)- approximated in 
O (m 4//3 e -3 ) time. 
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4.3 The Proof of Lemma [4.11 



All that remains is to prove the bounds given by Lemma 4.1 on the cardinality and capacity of H. 
To do so, we will use the effective resistance of the circuits on which we compute electrical flows 
as a potential function. The key insight is that we only cut an edge when its flow accounts for 
a nontrivial fraction of the energy of the electrical flow, and that cutting such an edge will cause 
a substantial change in the effective resistance. Combining this with a bound on how much the 
effective resistance can change during the execution of the algorithm will guarantee that we won't 
cut too many edges. 

Let r J be the resistances used in the j th electrical flow computed during the execution of the 
algorithrrj^J If an edge e is not in E or if e has been added to H by step j, set r{ = oo. We define 
the potential function 



<5>(j) = R eS ( r i)=£ r3 (f, 



where f r j is the (exact) electrical flow of value 1 arising from r J . Lemma 4.1 will follow easily from: 
Lemma 4.4. Suppose that F < F* < mF. Then: 

1. <5(j) never decreases during the execution of the algorithm. 

2. $(1) > m" 4 F- 2 . 

3. If we add an edge to H between steps j — 1 and j, then (1 — ^)^(j) > — !)• 



Proof. Proof of (1) 

The only way that the resistance ri of an edge e can change is if the weight w e is increased by 
the multiplicative weights routine, or if e is added to H so that ri is set to oo. As such, the 
resistances are nondecreasing during the execution of the algorithm. By Rayleigh Monotonicity 
(Corollary |2.5[ ), this implies that the effective resistance is nondecreasing as well. 

Proof of (2) 

In the first linear system, H = and r\ = 1- ^ 3 for all e E E. Let (5, V\S) be the minimum s-t cut 
of G. By the Max Flow-Min Cut Theorem (|1U| 15]). we know that the capacity u(S) = Y2eeE(S) Ue 
of this cut is equal to F* . In particular, 

! 1 + e/3 ^ 1 + e/3 1 

r = — > — > 

e u\ ~ F* 2 F* 2 

for all e € E(S). As f r i is an electrical s-t flow of value 1, it sends 1 unit of net flow across 
(5, V \ S); so, some edge e' G E(S) must have f r i(e') > 1/m. This gives 

$(1) = £AM = Y.^{e?rl > fAe'frl > (11) 
Since F* < mF by assumption, the desired inequality follows. 

Note that r J is not just the set of resistances arising from w 3 , since a single call to the oracle may compute 
multiple electrical flows as edges are added to H. 
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Proof of (3) 



Suppose we add the edge h to H between steps j — 1 and j. We will show that h accounts for a 
substantial fraction of the total energy of the electrical flow with respect to the resistances r-? , 
and our result will then follow from Lemma 12.61 

Let w be the weights used at step j — 1, and let / be the flow we computed in step j — 1. 
Because we added h to H, we know that congj(/i) > p. Since our algorithm did not return "fail" 

after computing this /, we must have that 



(12) 



Using the definition of 1 , the fact that congj(/i) > p, and Equation (12), we obtain: 

J{h?rl l = J{hf' 
>f[hV 



i l">li 

^2 e l w |l 
3mu 2 




to i 



> 



3(l + e)m 



congj(/i) 2 ((l + e)|to| 
£ rJ -i(7). 



3(1 + e)m j 



The above inequalities establish that edge h accounts for more than a 3 ^ e ^ m fraction of the 

total energy £ r i{f) of the flow /. 

The flow / is the approximate electrical flow computed by our algorithm, but our argument 
will require that an i nequ ality like this holds in the exact electrical flow f r j-i. This is guaranteed 



by part b of Theorem 
gives us that 



2.3 



rJ- 



which, along with the facts that £ r (f) > Srifri- 1 )^ P — x i an d e < 1/2, 

6/3 



2mR 

The result now follows from Lemma |2.6 



3(1 + e)m 



) ZrifrJ- 1 ) > Jzi^rifri-' 
2mK J om 



□ 



We are now ready to prove Lemma 4.1 



Proof of Lemma 4-1 Let fc be the cardinality of the set H at the end of the algorithm. Let / be 
the flow that was produced by our algorithm just before fc-th edge was added to H, let j be the 
time when this flow was output, and let to be the corresponding weights. 

As the energy required by an s-t flow scales with the square of the value of the flow, 



F 2 



F 2 



(13) 
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By the construction of our algorithm, it must have been the case that £ r j{f) < (1 + e)|to|i. This 
inequality together with equation ( 13 ) and part [2] of Lemma 4.4 implies that 



m = MM < £r jf [ < (i + e )l^|im 4 d>(i). 

Now, since up to time j we had k— 1 additions of edges to H, parts [T] and [3] of Lemma 4.4, and 
Lemma |3.4| imply that 



1 - — J < < (1 + e)|«;|im 4 < (1 + e)m 4 fmexp (ii±^jv)) < 2m 5 exp(3e _1 mm), 



5m / 3>(1) 

where the last inequality used the fact that e < 1/2. Rearranging the terms in the above inequality 
gives us that 

In 2 + 51nm + 3e _1 lnm 6e _1 lnm 30m In m 

In f 1 - %£) In(l-f-) e V 

V 5m y V 5m / 

where we used the inequalities e < 1/2 and log(l — c) < — c for all c S (0, 1). This establishes our 
bound on cardinality of the set H. 

To bound the value of u(H), let us note that we add an edge e to H only when we send at least 
pu e units of flow across it. But since we never flow more than F units across any single edge, we 
have that u e < F/p. Therefore, we may conclude that 

, r ^,F 30mFlnm 

u{H) < H - < — , 

p e^p 6 

as desired. □ 
4.4 Improving the Running Time to O (mn 1 ' 3 r 11 '' 3 ) 

We can now combine our algorithm with existing methods to further improve its running time. In 
|12j (see also [5]), Karger presented a technique, which he called "graph smoothing", that allows 
one to use random sampling to speed up an exact or (1 — e)-approximate flow algorithm. More 
precisely, his techniques yield the following theorem, which is implicit in [12 and stated in a more 
similar form in [5]: 

Theorem 4.5 ( [T2"l 15] ) . Let T(m,n,e) be the time needed to find a (1 — e)- approximately maximum 
flow in an undirected, capacitated graph with m edges and n vertices. Then one can obtain a 
(1 — e)- approximately maximal flow in such a graph in time O (e 2 m/n ■ T(0 (ne -2 ) , n, f2(e)) 

By applying the above theorem to our O (m 4 / 3 e~ 3 ) algorithm, we obtain our desired running 
time bound: 

Theorem 4.6. For any < e < 1/2, the maximum flow problem can be (1 — e)- approximated in 
O (mn'^e" 11 ' 3 ) time. 
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4.5 Approximating the Value of the Maximum s-t Flow in Time (m + n^ 3 e 8 / 3 ) 

Given any weighted undirected graph G = (V, E, w) with n vertices and m edges, Benczur and 
Karger [3] showed that one can construct a graph G' = (V, E', w') (called a sparsifier of G) on the 
same vertex set in time O (m) such that \E'\ = 0(n log n/e 2 ) and the capacity of any cut in G' is 
between 1 and (1 + e) times its capacity in G. Applying our algorithm from Section [4] to a sparsifier 
of G gives us an algorithm for (1 — e)-approximating the value of the maximum s-t flow on G in 
time O (m + n 4 / 3 e -3 ) . 

We note that this only allows us to approximate the value of the maximum s-t flow on G. It 
gives us a flow on G' , not one on G. We do not know how to use an approximately maximum s-t 
flow on G' to obtain one on G in less time than would be required to compute a maximum flow in 
G from scratch using the algorithm from Section |4j 

For this reason, there is a gap between the time we require to find a maximum flow and the time 
we require to compute its value. We note, however, that this gap will not exist for the minimum 
s-t cut problem, since an approximately minimum s-t cut on G' will also be an approximately 
minimum s-t cut G. We will present an algorithm for finding such a cut in the next section. By the 
Max Flow-Min Cut Theorem, this will provide us with an alternate algorithm for approximating 
the value of the maximum s-t flow. It will have a slightly better dependence on e, which will allow 
us to approximate the value of the maximum s-t flow in time O (m + n 4//3 e -8 / 3 ). 



5 A Dual Algorithm for Finding an Approximately Minimum s-t 
Cut in Time O (m + n 4 / 3 e~ 8 / 3 ) 

In this section, we'll describe a dual perspective that yields to an even simpler algorithm for 
computing an approximately minimum s-t cut. Rather than using electrical flows to obtain a flow, 
it will use the electrical potentials to obtain a cut. 

The algorithm will eschew the oracle abstraction and multiplicative weights machinery. Instead, 
it will just repeatedly compute an electrical flow, increase the resistances of edges according to the 
amount flowing over them, and repeat. It will then use the electrical potentials of the last flow 
computed to find a cut by picking a cutoff and splitting the vertices according to whether their 
potentials are above or below the cutoff. 

The algorithm is shown in Figure [6| It finds a (1 + e)-approximately minimum s-t cut in time 
O (m 4//3 e _8//3 ); applying it to a sparsifier will give us: 

Theorem 5.1. For any < e < 1/7, we can find a (1 + e) -approximate minimum s-t cut in 
O (m + n 4 / 3 e~ 8 / 3 ) time. 

We note that, in this algorithm, there is no need to deal explicitly with edges flowing more than 
p, maintain a set of forbidden edges, or average the flows from different steps. We will separately 
study edges with very large flow in our analysis, but the algorithm itself avoids the complexities 
that appeared in the improved flow algorithm described in Section |4j 

We further note that the update rule is slightly modified from the one that appeared earlier in 
the paper. This is done to guarantee that the effective resistance increases substantially when some 
edge flows more than p, without having to explicitly cut it. Our previous rule allowed the weight 
(but not resistance) of an edge to constitute a very small fraction of the total weight; in this case, 
a significant multiplicative increase in the weight of an edge may not produce a substantial change 
in the effective resistance of the graph. 
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Input : A graph G = (V,E) with capacities {u e } e , and a target flow value F 
Output: A cut (S,V\S) 

Initialize w° <- 1 for all edges e, p <- 3m 1 / 3 e -2 / 3 , N <- he^^m 1 ^ km, and 5 <- e 2 
for j := 1, . . . , N do 

Find an approximate electrical flow /* and potentials using Theorem 



2.3 



on G with 



resistances 1 = " a , target flow value F, and parameter (5. 

Wg <— Wg _1 + |cong^_! {e)w l ~ 1 + 1 f° r eacn e £ E 



Scale and translate so that S , = 1 and cj> t — 
Let 5a, = € V | > a;} 

Set S to be the set S x that minimizes (S^, V \ S x ) 

If the capacity of (S x , V \ S x ) is less than F/(l - 7e), return (S X) V \ S x ). 
end 

return "fail" " 



Figure 6: A dual algorithm for finding an s-t cut 



5.1 An Overview of the Analysis 

To analyze this algorithm, we will track the total weight placed on the edges crossing some minimum 
cut. The basic observation for our analysis is that the same amount of net flow must be sent across 
every cut, so edges in small cuts will tend to have higher congestion than edges in large cuts. 
Since our algorithm increases the weight of an edge according to its congestion, this will cause our 
algorithm to concentrate a larger and larger fraction of the total weight on the small cuts of the 
graph. This will continue until almost all of the weight is concentrated on approximately minimum 
cuts. 

Of course, the edges crossing a minimum cut will also cross many other (likely much larger) 
cuts, so we certainly can't hope to obtain a graph in which no edge crossing a large cut has non- 
negligible weight. In order to formalize the above argument, we will thus need some good way to 
measure the extent to which the weight is "concentrated on approximately minimum cuts" . 



In Section 5.2 we will show how to use effective resistance to formulate such a notion. In 
particular, we will show that if we can make the effective resistance large enough then we can find 
a cut of small capacity. In Section |5.3| we will use an argument like the one described above to 
show that the resistances produced by the algorithm in Figure ^ converge after N = O (m 1 / 3 e -8 / 3 ) 
steps to one that meets such a bound. 



5.2 Cuts, Electrical Potentials, and Effective Resistance 

During the algorithm, we scale and translate the potentials of the approximate electrical flow so 
that (p s = 1 and <f> t = 0. We then produce a cut by choosing x G [0, 1] and dividing the graph into 
the sets S = {v £ V \ tp v > x} and V \ S = {v £ V \ (f) v < x} . The following lemma upper bounds 
the capacity of the resulting cut in terms of the electrical potentials and edge capacities. 

Lemma 5.2. Let <fi be as above. Then there is a cut S x of capacity at most 

X] \<f>u-<i>v\ u (u,v)- ( 14 ) 
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Proof. Consider choosing x E [0, 1] uniformly at random. The probability that an edge (u, v) is cut 



is precisely \<f>(u) — <f>(v)\. So, the expected capacity of the edges in a random cut is given by (14) 



and so there is a cut of capacity at most ( 14 ) . □ 



Now, suppose that one has a fixed total amount of resistance fj, to distribute over the edges of 
a cut of size F. It is not difficult to see that the maximum possible effective resistance between s 
and t in such a case is and that this is achieved when one puts a resistance of y on each of the 



edges. This suggests the following lemma, which bounds the quantity in Lemma 5.2 in terms of the 
effective resistance and the total resistance (appropriately weighted when the edges have non-unit 
capacities): 

Lemma 5.3. Let fi = Y2e u e r e> an d ^ ^ e effective s-t resistance of G with edge resistances given 
by r be R c ff(r). Let <f> be the potentials of the electrical s-t flow, scaled to have potential drop 1 
between s and t. Then 



v 'Mr) 



If, (f) is an approximate electrical potential returned by the algorithm of Theorem 2.3 when run with 
parameter 5 < 1 /3, re-scaled to have potential difference 1 between s and t, then 



g* eK <(l + M)^. 



Proof. By Fact 2.4 the rescaled true electrical potentials correspond to a flow of value l/i? e ff( r ) 
and 

^ r e R eS {r)' 
So, we can apply the Cauchy-Schwarz inequality to prove 



ulr, 



e' e 



By parts a and c of Theorem 2.3 after we rescale </> to have potential drop 1 between s and t, it 
will have energy 

v 0(e) 2 1 + 8 1 1 
4" r e ~ l-8R cS (r) ~ { ^ 1 i? cff (r)' 

for 5 < 1/3. The rest of the analysis follows from another application of Cauchy-Schwarz. □ 

5.3 The Proof that the Dual Algorithm Finds an Approximately Minimum Cut 

We'll show that if F > F* then within N = be~ & ^m 1 ^ In to iterations, the algorithm in Figure [H] 
will produce a set of resistances r % such that 



(Ff 



Resir^y (1-76)^-2. (15) 
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Once such a set of resistances has been obtained, Lemmas 5.2 and 5.3 tell us that the best potential 
cut of <j> w iU have capacity at most 

1 + 25 ^ 1 „ 
--F < —F. 



VT^Te ~ 1 - 7e 

The algorithm will then return this cut. 

Let C be the set of edges crossing some minimum cut in our graph. Let uc = F* denote the 
capacity of the edges in C. We will keep track of two quantities: the weighted geometric mean of 
the weights of the edges in C, 

(\ l/«c 
eeC / 



and the total weight 



of the edges of the entire graph. Clearly v l < max e ^c w l- I n particular, 

v* < H* 

for all i. 

Our proof that the effective resistance cannot remain large for too many iterations will be 
similar to our analysis of the flow algorithm in Section |4j We suppose to the contrary that Reff * t < 
(1 — 7e)jj^2 f° r each 1 < i < N. We will show that, under this assumption: 



1. The total weight // doesn't get too large over the course of the algorithm [Lemma 



5.4 



2. The quantity v % increases significantly in any iteration in which no edge has congestion more 



than p [Lemma 5.5 . Since p} doesn't get too large, and u l < p\ this will not happen too 
many times. 

3. The effective resistance increases significantly in any iteration in which some edge has conges- 



tion more than p [Lemma 5.6 . Since p l does not get too large, and the effective resistance 



is assumed to be bounded in terms of the total weight p l , this cannot happen too many times. 
The combined bounds from [2] and [3] will be less than N, which will yield a contradiction. 
Lemma 5.4. For each i < N such that R c g(r l ) < (1 — 7e)-^2, 

P < P exp ' 



P 

So, if for all i < N we have R c ff(r l ) < (1 — 7e)-^, then 



p N <p°e W ( €{1 ~ 2e) N). (16) 
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Proof. If R e g(r l ) < (1 — 7e)j^ then the electrical flow / of value F has energy 

F 2 ^ cfT (r i ) = ^con g/i (e) 2 «4 < (1 - 7e)/A 



By Theorem 2.3, the approximate electrical flow / has energy at most (1 + 5) times the energy of 
/, 

]Tcon g/ ~( e ) 2 < < (l + <5)(l-7e)^ < (1 - 6e)/A 
Applying the Cauchy-Schwarz inequality, we find 



]Tcon g/ ~(eK < ^^^con g/ ~(e)2^ < vT=6i/i* < (1 - 3e)^ 
Now we can compute 



e 2 



^Z W e + -X^ COn g/'( e ) + M 



<|1 + ^'|M' 



< // exp 



e(l-2e) 



as desired. □ 
Lemma 5.5. If congp (e) < p for all e, then 

S» > exp ( £(1 - £) ) „\ 

Proof. To bound the increase of z^ +1 over z/, we use the inequality 

(1 + ex) > exp (e(l — e)x) , 
which holds for e and x between and 1. We apply this inequality with x = con g j t (e)/p. As /' is 
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flow of value F and C is a cut, X^eeC fl 



> F. We now compute 

l/uc 



VeeC / 

=^ (n( i +^( e )) 



> z/ 



ex p ( ~~ X] 



u e — - — cong^(e) 



1 ^e(l-e) 



> v l exp 

> v % exp 



1 e(l - e) 

e(l-e) 
P 



□ 



Lemma 5.6. If R e s(r l ) < (1 — 7e)j^ and there exists some edge e such that congr^e) > p, t/ien 



R c s(r i+1 ) > exp 
Proof. We first show by induction that 



4m 



If i = we have 1 > e. For i > 0, we have 



>-(- + -) 



e /e(l-2eV 

> — exp j /i- 

m \ P 

> ~p i+ \ 
m 



by Lemma 5.4 
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We now show that an edge e for which congj 4 (e) > p contributes a large fraction of the energy 
to the true electrical flow of value F. By the assumptions of the lemma, the energy of the true 
electrical flow of value F is 

F 2 R eS {r) < (1 - 7e)/A 
On the other hand, the energy of the edge e in the approximate electrical flow is 

f{e) 2 ri = cong ? (e) 2 w e > p 2 —p\ 

As a fraction of the energy of the true electrical flow, this is at least 

1 p 2 e _ 1 
l-7e m ~ (1 - 7e)e 1 /3 m i/3- 



By part 6 of Theorem 2.3 the fraction of the energy that e accounts for in the true flow is at least 

1 e*_ 1 _ ph 

(1 - 7e)e 1 /3 m i/3 2mR ~ e^m 1 / 3 ~ m ' 

As 

w\ +1 > w\ (l + ^cong^e)^ > (1 + e)u& 



we have increased the weight of edge e by a factor of at least (1 + e). So, by Lemma 2.6 



□ 

We now combine these lemmas to obtain our main bound: 

Lemma 5.7. For e < 1/7, after N iterations, the algorithm in Figure^ will produce a set of 
resistances such that R e g(r l ) < (1 — 7e) j^. 



Before proving Lemma |5.7[ we note that combining it with Lemmas 5.2 and 5.3 immediately 
implies our main result: 

Theorem 5.8. On input e < 1/7, the algorithm in Figure^runs in time O (m 4 / 3 e~ 8 / 3 ) . If F > F* , 
then it returns a cut of capacity at most F/(l — 7e), where F* is the minimum capacity of an s-t 
cut. 

To use this algorithm to find a cut of approximately minimum capacity, one should begin as 
described in Section [2 .2| by crudely approximating the minimum cut, and then applying the above 
algorithm in a binary search. Crudely analyzed, this incurs a multiplicative overhead of 0(log m/e). 



5.4 



Proof of Lemma \5. 7\ Suppose as before that R e g (r l ) < (1 — 7e) for all 1 < i < N. By Lemma 
and the fact that pP = m, the total weight at the end of the algorithm is bounded by 

< p° exp ( 6(1 ~ 2e M = mexp ( ^ ~ 2e) N ) . (17) 
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Let A C {1, . . . , N} be the set of i for which cong^(e) < p for all e, and let B C {1, . . . , iV} 
be the set of i for which there exists an edge e with congj^e) > p. Let a = \ A\ and b = \B\, and 
note that a + b = N . We will obtain a contradiction by proving bounds on a and 6 that add up to 
something less than N. 

We begin by bounding a. By applying Lemma 5.5 to all of the steps in A and noting that v l 
never decreases during the other steps, we obtain 



v > exp 



e(l-e) \ o 

a I v = exp 



el 



(18) 



Since i/ < /u , we can combine Equations (17) and (18) to obtain 

/'e(l-e) \ /e(l-2e) Ar 
exp I — -aj < mexp I — -N 

Taking logs of both sides and solving for a yields 

e (l_ € ) € {l-2e) 



a < 



P 



N + lnm, 



from which we obtain 
1 - 2e 



a < 



1 



N + 



6(1 - e) 



lnm < (1 - e)N + 



7p 



el 



— lnm < (1 — e)A^ H In m. 

e) 6e 



(19) 



We now use a similar argument to bound b. By applying Lemma 5.6 to all of the steps in B 



and noting that R c s(r % ) never decreases during the other steps, we obtain 



R e s(r N ) > exp 



4m 



b R cS (r°) > exp 



eV 
4m 



1 + e 

2 ^?*2 



??? 



> exp 



eV 2 
-. — i 
4m 



1 + e 



m 



where the second inequality applies the bound i? e g(r°) > l/(m 2 -F* 2 ), which follows from an argu- 
ment like the one used in Equation (11). Using our assumption that R e g(r ) < (1 — 7e)^2 an d 
the bound on p, N from Equation (17), this gives us 



eV, 
exp | — — o 



1 + e 



< 



1 - 7e 



4m J m 2 F 2 F 2 
Taking logs and rearranging terms gives us 

1 -7e 



mexp 



e(l - 2e) 



N 



4m(l - 2e) 4m , 
b < i— = -iV + In 



ep» 



e 2p2 



q\ 4m , T 12m , 

ttT < -^-iV + lnm. 

1 + e J ep 6 e z p A 



(20) 



Adding the inequalities in Equations (19) and (20), grouping terms, and plugging in the values of 
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p and N, yields 

f,-, n 4m \ t» t /7p 12m \ 
^ = a + & <^(l- e) + ^j7V + ^ + _jw 

= (1 - e) + — ^ (5e~ 8/3 m 1/3 lnm) + h ^7 ^ — 77777 lnm 

V V ; e(27me~ 2 ) J y ' I 2e e 2 (9m 2 / 3 e" 4 /3) J 

n 4 A / r -8/3 1/3! x /Tm 1 ^ ^mV'A 
= (^(1 - 6) + Y7 ) (5e 8 / 3 mV3l nm) + ^_ + __j l nm 

5 m l/3 / 41 12e \ m l/3 



This is our desired contradiction, which completes the proof. □ 
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A Computing Electrical Flows 



Proof of Theorem \2.3\ Without loss of generality, we may scale the resistances so that they lie 
between 1 and R. This gives the following relation between F and the energy of the electrical s-t 
flow of value F: 

F 2 

— < £ r (f) < F 2 Rm. (21) 
m 

Let L be the Laplacian matrix for this electrical flow problem. Note that all off-diagonal entries 
of L lie between —1 and —1/R. Recall that the vector of optimal vertex potentials <f> of the electrical 
flow is the solution to the following linear system: 



L(f> = F Xs ,t 

For any e > 0, the algorithm of Koutis, Miller and Peng [13] provides in time O (m log 2 m log 2 log m log e _1 ) 
a set of potentials such that 

ii - " < £ 



where 



Note also that 



4,-4, 

U\\ L 



L 

def 



L i 



cf) T Lcf). 
£r(f), 
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where / is the potential flow corresponding to 0. Now, let / be the potential flow corresponding 
to 4>, 

f = CB T $. 

The flow / is not necessarily an s-t flow because only approximately satisfies the linear system. 
A small amount of flow can enter or leave every other vertex. In other words,^f may not satisfy 
Bf = Fxs,t- Below, we will compute an approximation / of / that satisfies Bf = F\s,t- 

Before fixing this problem, we first observe that the energy of / is not much more than the 
energy of / . We have 







Srif), 



and by the triangle inequality 







<II0IU + 



0-0 



<(l + e)U0ll 



L ■ 



So 



£r(f)<(l + e) 2 £ r (f). 



To see that the flow / does not flow too much into or out of any vertex other than s and t, note 
that the flows into and out of vertices are given by the vector 

iext = Bf. 

Note that the sum of the entries in i ex t is zero. We will produce an s-t flow of value F by adding 
a flow to / to obtain a flow / for which 

Bf = F X s,t. 

Let 77 be the maximum discrepancy between i ext and Fxs,t- 

def I, . _ „ || 
V — II "text ~ " Xs,t II oo • 



We have 

V < \\iext ~ FXa,t\\ 2 



L<t>-L<t> 



< 11X11 



0-0 



< 2n 



0-0 



< 2ne v / £ r (f). 



It is an easy matter to produce an s-t flow / that differs from / by at most nr/ on each edge. 
For example, one can do this by solving a flow problem in a spanning tree of the graph. Let 
T be any spanning tree of the graph G. We now construct a flow in T that with the demands 
FXs,t{u) — i e xt{u) at each vertex u. As the sum of the positive demands in this flow is at most nr], 
one can find such a flow in which every edge flows at most nrj. Moreover, since the edges of T are 
a tree, one can find the flow in linear time. We obtain the flow / by adding this flow to /. The 
resulting flow is an s-t flow of value F. Moreover, 



/-/ 



< nrj. 



To ensure that we get a good approximation of the electrical flow, we now impose the requirment 
that 

5 

6 - 2n 2 mV2 jR 3/2 • 
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To check that requirement a is satisfied, observe that 

e 

< J2 r e{fe + nri) 2 

e 

< £ r (f) + 2nr, r e / e + n 2 rf r e 

e e 

< £ r (f) + {2nr)F + n 2 7] 2 ) r e 

e 

< £ r {f) + e(2n 2 {mR) 1/2 + An A )£ r (f)mR 
<£ r {f) ((1 + e) 2 + e6n 4 mi? 3/2 ' 



We may ensure that this is at most (1 + 5)£ r (f) by setting 

6 



To establish part b, note that 



12n 4 mR 3 / 2 ' 



and so 



We now bound 



4>~4> 



< e 2 £ r (f). 



Tefi-TefA by 



r e f e - r ef e = r e {fe ~ fej r e (/ e + / e 

<e 2 £ r (f)(2 + e) 2 £ r (f), 



which implies 



It remains to bound 



r e fe -refl < e(2 + e)£ r (f) . 
, which we do by the calculation 



fe fe 



fe + fe 



r efe r ef e 

< v / i?2n 2 e v /f r (/)(2 + e)y/£Afj 
= e\ r R2(2 + e)n 2 £ r (f). 

Putting these inequalities together we find 



r ef e r ef e 



(by 



< eVR(2n 2 + 1)(2 + e)£ r (f) < £ r {f) 



2mR 



To prove part c, first recall that 



L<j> = F Xs 
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£ r (f) = F 2 R cS (r). 



e 2 H\\l> 



L 



(4-$) l(4>-$) 



<t> 



+ 



2 ~ T 

L -2(f> Lcf) 



\\<t>\\l + 



<t> 



2Fxl t 4>>\\<t>\\ 2 L + 



4> 



\\4>\\] 



>(l + (l-6) 2 -e 2 ) ||*||] 

= (2-2 e )||^||| 

= (2-2e) F 2 R cS (r). 
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